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Abstract: We study a supersymmetric extension of the vector-like lepton scenario, such 
that the vacuum instability induced by large lepton Yukawa couplings is lifted by the 
presence of superpartners at or below the TeV scale. In order to preserve the unification 
of gauge couplings, we introduce a full 16 + 16 of SO(10), and determine the maximal 
possible values for the Yukawa couplings consistent with perturbativity at the GUT scale. 
We find that the Higgs to diphoton decay rate can be enhanced by up to 50% while 
maintaining vacuum stability and keeping the new particle masses above 100 GeV, while 
larger enhancements are possible if the masses of the new particles are lowered further. 



1 Introduction 

In Ref. [1] we have presented a model where the Higgs decay rate to diphotons is increased 
through loops of mixed vector-like leptons. A vector-like doublet and a vector-like singlet 
allow for both Yukawa and Dirac masses. The resulting mixing leads to a sign flip of the 
coupling of the lightest lepton to the Higgs, such that constructive interference with the 
standard model (SM) amplitude for h — > 77 is possible, resulting in a diphoton decay rate 
that is enhanced relative to the SM. Similar models were presented in [2, 3], and further 
studies have appeared in [4-17] . 

One of the shortcomings of the model [1] is that the additional 0(1) Yukawa couplings 
accelerate the downward running of the Higgs quartic coupling, such that vacuum stability 
is lost at scales around 10 TeV. This tension between an enhanced di-photon rate and 
vacuum stability was first noted in [1, 2] and further explored in [23-28]. The goal of the 
present paper is to provide a UV completion of the model [1], which guarantees vacuum 
stability up to very high scales and unification of gauge couplings at roughly 10 16 GeV. 

In supersymmetric models, the scalar potential is fixed by gauge invariance and su- 
persymmetry relations, so the RGE evolution of the quartic coupling does not lead to a 
stability problem, provided that it is well behaved at scales below the SUSY breaking scale. 
It is therefore natural to try to embed the model [1] into a supersymmetric extension of 
the SM. 

It is well known that within the MSSM the gauge couplings unify at approximately 
10 16 GeV with very high accuracy. In order to not destroy the sensible relation between the 
values of the couplings at the low scale and the beta function coefficients, additional matter 
fields have to be embedded in complete multiplets of SU(5). A vector-like doublet with 
leptonic quantum numbers is naturally contained in a 5 + 5 representation of SU(5), while 
a corresponding singlet field can emerge from a 10 + 10. In the language of SO (10) grand 
unification, these representations can be combined into a 16 + 16 representation where the 
additional degrees of freedom have neutrino-like quantum numbers. 

Additional matter with couplings to the Higgs sector will also give a positive contribu- 
tion to the mass of the lightest Higgs boson in the MSSM [29-33], thus might help reducing 
the fine tuning problem that is present in the Higgs sector of the MSSM. 

This paper is organized as follows: In Sec. 2, we introduce the model and the particle 
content at the weak scale. The RGE evolution of gauge and Yukawa couplings are discussed 
in Sec. 3. The effects of the new particles on the Higgs mass and di-photon decay rate are 
calculated in section 4, while in Sec. 5 the conditions for vacuum stability are derived. In 
Sec. 6 we present numerical results before concluding in Sec. 7. 

2 The Model 

Our model is the minimal supersymmetric standard model (MSSM) extended by a 16 + 16 
of SO(10), such that the unification of gauge couplings at the GUT scale is guaranteed. 
Models with such a particle spectrum were previously studied e.g. in [29] [30]. Below Mqut 

1 Effects of vector-like quark multiplets are discussed e.g. in [18-22]. 
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E'L 


AT- 


Quantum Numbers 


(2, -1/2) 


(1, 1) 


(1, o) 


(2, 1/2) 


(1, -1) 


(1,0) 



Table 1. Additional fields below the TeV scale. All fields are left-handed superfields, and the 
quantum numbers specify the transformations of the fields under the SU(2)xU(l) gauge group of 
the SM. Primes indicate fields coming from the 16 multiplet, while the double primed fields originate 
from the 16 multiplet. 



the SO(10) multiplets are split. We assume that all the additional colored states obtain 
masses above the TeV scale, and therefore are in agreement with current LHC limits. 

Below the TeV scale, we assume that the only fields beyond the MSSM particle content 
are the un-colored components of the original 16 + 16 supermultiplets. In terms of chiral 
superfields, these are SU(2) doublets L' L and L R , the singlet fields E' R and E L as well as 
singlet neutrino superfields N' R and N L with the quantum numbers indicated in Tab. 1. 
Our notation is such that the bar extends over the implicit chiral projector, i.e. E' R is a 
left-handed superfield. The superpotential is 



W =W MS sm - M L L' L L" R + M E E'[E' R - y' c L' L H d E' R + y" c L" R H u E" L 



-y' n L' L H u N' R + y^L" R H d N'l 



MijNiNj , 



(2.1) 



where we have neglected the colored fields beyond the MSSM, and (Aq,^) = (N' R ,N L ) 
since the neutrinos can have Majorana mass terms in addition to the Dirac mass terms 
for L and E fields. As in [1], we impose a parity symmetry under which the vector-like 
multiplets are odd, such that mixing with MSSM leptons and sleptons is forbidden. 

Contraction of SU(2) indices is implicit in (2.1). In analogy with the leptonic su- 
perfields, the Higgs superfields are defined as H u = (H+,H°) T and H d = {H%HJ) T . 
Indices are contracted using the anti-symmetric tensor e with e\2 = 1, for example L' L H d = 
(N' L H^ — E' L H^). The sign conventions in (2.1) are chosen such that all entries in the 
charged fermion mass matrix come with a positive sign and all entries in the neutral 
fermion mass matrix with a negative sign, for convenience. 

Explicitly, the charged fermion mass matrix is given by 



e L e L e R e R 



Mi 



R 



+ h.c. with. Me 



VkJ 



( y' c v d M L \ 
M E y" c v u 
y' c v d M E 
\M L y" c v u / 



(2.2) 



where we used lower case letters for the fermionic components of the superfields, and 
we have introduced the scalar vacuum expectation values v u = (H®) and v d = (H®). 
In the supersymmetric limit, the corresponding slepton mass matrix is simply given by 
M 2 ^ = M e Me- The [xH u H d term communicates SUSY breaking to the lepton sector, 
such that, in a non-trivial Higgs background, the charged slepton mass matrix assumes the 
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form 



M\ = (2.3) 

(\y' c \ 2 v 2 + \M L \ 2 + m 2 e , L y'*v d M E + M L y" c v u -y' c v u »* b L 

M* E y' c v d + y't*v u M L \y^\ 2 v 2 u + \M E \ 2 + m 2 e „ b E -y"v d ^ 

-y'*v ul i b E \M E \ 2 + \y' c \ 2 v 2 + m 2 e , R y'*v d M L + M* E y^v u 

\ b L -y?v df i M* L y' c v d + y^v u M E \M L \ 2 + \y^\ 2 v 2 u + m 2 J 

in the basis (e' L , e L , e' R , e R ). The D-term contributions are not explicitly written down 
since their effects are small, but they are included in our analysis. Note that, in addition 
to the supersymmetric mass terms, we have allowed for the following bi-linear soft breaking 
terms in the potential: 



V soit/ = m\ \£' L \ 2 + m\ \i" R \ 2 + m 2 e ,,\e L \ 2 + m 2 Je' R \ 2 + \{b L l'lt" R + b E ~e'* R e" L + h.c.) . (2.4) 



The structure of the mass matrices and their effects on the h — > 77 rate will be further 
explored in Sec. 4. It is worth noting that half of the scalar degrees of freedom can be 
decoupled by increasing in rt and m „ to high values. This limit is similar to the light stau 
scenario [34, 35], with a 2 x 2 charged slepton mass matrix 

2 f\y' c \ 2 v 2 + \M L \ 2 + m 2 e , L -y' c v u ^ \ 

Z' 2 * 2 { -y'*v u » \M E \ 2 + \y' c \ 2 v 2 + m 2 e J- ^ 

A similar contribution is obtained when instead the double primed fields e'[ and e R are 
lifted, however in that case the mixing is proportional to v d and therefore suppressed in 
the large tan f3 regime. 

Before moving to the next section, let us briefly review existing experimental limits 
on uncolored charged scalars and fermions. The LEP experiments have searched for such 
particles, and their results are collected in the particle data booklet (PDG) [36]. 

Limits on additional charged leptons are given explicitly in the PDG. For a charged 
lepton that decays to W^is, a limit of 100.8 GeV is quoted, while a limit of 101.9 GeV is 
given for a charged lepton that is not degenerate with it's corresponding neutrino. However, 
as discussed in detail in [1], if a neutrino- like state is close in mass, these limits become 
invalidated. 

Limits on the scalar partners of standard model leptons strongly depend on the flavor 
composition of the sleptons, and range from 107 GeV for left-handed selectrons down to 
81.9 GeV for staus. As in the fermionic case, most of these bounds are weakened if other 
neutral states are close by in mass, and none of them can be directly applied to scalar 
partners of new leptons, as in our model. As absolute lower bound on the mass of the 
lightest leptons and slepton we therefore impose m E > m^/2 « 62.5 GeV. To facilitate 
the discussion in the remainder of the paper, we define two LEP limits: 

• Conservative LEP limit: m El > 100 GeV, m_g > 90 GeV, 
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Figure 1. Evolution of gauge couplings in the MSSM (dashed lines) and in the MSSM extended 
by a 16 + 16 of SO(10) (solid lines). From top to bottom the curves correspond to ai, a-i and a^. 

• Optimistic LEP limit: m El > 62.5 GeV, > 62.5 GeV, 

with the understanding that a spectrum that satisfies the conservative LEP limit will 
certainly be in agreement with current experimental bounds, while a spectrum that satisfies 
the optimistic LEP limit might require additional invisible neutral states to be close in mass 
in order to not be in conflict with existing searches. It should be noted that such states 
are present in our model in the form of new neutrinos and sneutrinos. 

The possibility for the LHC experiments to discover or constrain such particles is 
discussed for example in [2, 35], with the general consensus that high luminosity at the 
14 TeV LHC is required. 

3 Running of Couplings 

The renormalization group evolution (RGE) of the gauge and Yukawa couplings is strongly 
affected by the presence of additional matter charged under the strong and weak interac- 
tions. To simplify the analysis of this section, we assume a common threshold at the TeV 
scale for the new vector-like states and for all SUSY partners. While ultimately we will 
be interested in scenarios where some of the vector-like matter is lighter, the effects on the 
RGE of the gauge couplings are negligible. 

The one-loop evolution of gauge couplings is governed by 



ar\A) = a r\M z ) 




log 



A 



(3-1) 
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When A > 1 TeV, the one loop beta function coefficients are given by 

' I ™> ?• 

/3 2 = 5 - 2N f = -5 , 
fo = 9-2N f = -l, 

with Nf = 5, while the corresponding values in the MSSM are obtained by taking Nf = 3. 
Note in particular that this leads to a sign change in the beta function for a%, which 
shows that the full particle content of the MSSM +16 + 16 is enough to render the strong 
interactions asymptotically non-free. 

Fig. 1 shows the evolution up to high scales. Unification of gauge couplings occurs 
roughly around Mjj ~ 1.5 x 10 16 GeV, with coupling strengths of about an ~ 0.15 corre- 
sponding to gi ~ 1.4. At the one loop level this is compatible with perturbative unification. 
For a discussion of higher order effects see e.g. [30]. 

Yukawa couplings have a well known fixed point behavior in the infrared, which allows 
us to determine an upper limit on the magnitude of the couplings at the electroweak scale. 
At the one-loop level, the RGE equations for the Yukawa couplings take the form 

= ~Mm) . ( 3 - 2 ) 
an n 

with beta functions given by [37] 

MM) = ^jW (6y t 2 + yl - 4vr a x + 3a 2 + f a 3 )) , (3.3) 
1 

167r 2; 

& 2 (M) = ^2/2 (4yl + 3y t 2 - 4vr (fa! + 3a 2 )) , (3.5) 

where we have suppressed the scale dependence on the right-hand sides. For sufficiently 
large initial values at the weak scale, yi and y 2 will diverge at high scales. The top Yukawa 
coupling is yt(Mz) = Mt/(vs'mf3), such that the bound on y 2 will be more stringent for 
small tan/3, while for larger tan/3 the bottom and tau Yukawas will lead to a slightly 
stronger bound on y\. 

For the numerical analysis, we used the MS running masses of the top and bottom 
quarks at the weak scale, M t {M z ) ~ 165 GeV and M b (M z ) ~ 2.7 GeV. Then the top and 
bottom Yukawas are given by 

Y,m)= M ^f, (3.6, 
v smp 

n(Mz) - vco S (3 1 + AilV (3 - 7) 

The SUSY-QCD correction to the bottom mass, AMb, is relevant in the large tan/3 regime 
and was included in our analysis. 



W = TTT^ 1 ( 4 ^i " 47r (l ai + 3a2 )) ' ( 3 - 4 ) 
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Figure 2. RGE evolution of Yukawa couplings for tan/3 = 2 (left) and tan/3 = 60 (right). Shown 
are the evolution of the new lepton Yukawa couplings Y' and Y" for various input values at the 
high scales, to illustrate the fixed points. 



The running of the couplings for small and large tan /3 is shown in Fig. 2. At the 
weak scale, the couplings y%, yi are largely independent of their precise magnitude at the 
unification scale. In particular, we note that Yukawas of order 0.5-0.9 at the weak scale are 
natural given order one input values at the high scale. For tan (3 ~ 60 the bottom Yukawa 
coupling is comparable to Y t , and the upper bounds on the new Yukawas are y' c < 0.9 and 
He ~ 0-8' while for small values of tan f3 ~ 2 the upper bounds are approximately given by 
y' c < 0.94 and < 0.72. 

Comparing with the running of the top and bottom Yukawa coupling, we note that 
the as contribution tends to suppress the quark Yukawas at high scales, or, reversing the 
argument, lead to larger couplings at low scales. If one would like to obtain larger y\^ at 
the weak scale, one could therefore either give up on perturbativity up to the unification 
scale or introduce additional interactions at intermediate scales that enhance the leptonic 
Yukawa couplings. 

Before moving to the next section, let us briefly comment on the evolution of the Higgs 
quartic couplings below the SUSY breaking scale. At the weak scale, the particle content 
of our model is that of the standard model with an additional Higgs doublet, a set of 
vector-like leptons as in [1], and potentially their super partners. We have seen in [1] that 
additional order one Yukawa couplings in conjunction with a Higgs mass of 125 GeV lead 
to vacuum instabilities due to the RGE evolution of the Higgs quartic coupling, but only 
above the TeV scale. Supersymmetry is expected to be restored at least partially around 
that scale 2 , above which the quartic couplings in the Higgs potential are determined trough 
supersymmetric relations and stop being a threat to the stability of the vacuum. 



4 Higgs Properties 

The MSSM Higgs sector consists of the two Higgs doublets H u , that acquire vacuum 
expectation values (VEV) v u and Vd, respectively. The VEVs are subject to the constraint 

2 In particular stops and weak gauginos are still allowed to be significantly lighter than one TeV. 
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v u + v d = (^74: GeV) 2 and are therefore parametrized as v u = v sin/3 and v d = v cos /3. 

The physical spectrum contains two neutral CP even Higgs bosons h° and H°, a 
neutral CP odd boson A and a charged scalar . The neutral mass eigenstates are linear 
combinations of the neutral doublet components H®, H®, with the mixing parameterized 
as: 

H u\ _ ( v u\ J_( cosa sincA f h°\ i_ ( cos p sin/3 \ /g°\ 
H° d ) ~ \y d ) + 71 \-sma cosa) [^H J + 71 l-sin/3 cos/3 J \A° ) ' 1 '"' 

Here a is the CP-even mixing angle, G° is the neutral Goldstone boson that can be removed 
by going to unitary gauge, and by convention h is defined to be the lightest of the CP-even 
mass eigenstates. 

A particularly interesting regime is the so called decoupling limit, which is characterized 
by a large mass scale txia for the non-standard Higgs bosons H°, A and . In this limit, 
the lightest CP even Higgs h° becomes SM like, and the mixing angle a is determined by 
a = (3 — 7r/2. Supersymmetric constraints on the potential constrain its mass to be 

m 2 = m% cos 2 (2/3) + radiative corrections. (4-2) 

For the reminder of this paper, we will work in the decoupling limit with M4 = 1 TeV, 
and further assume that the MSSM superpartners do not significantly modify the Higgs 
properties. This is strongly motivated by the recent discovery of a scalar resonance with 
a mass of around 125 GeV and with properties consistent with a SM Higgs boson. While 
this mass exceeds the lightest Higgs mass in the MSSM at the tree level, it is well known 
that radiative corrections can bring this mass in agreement with the observation. 

The dominant hP production channel at the LHC is gg — > h°, where the maximum 
contribution comes from the top loop. As new leptons and sleptons added to the MSSM 
are not colored, the production channel is not affected significantly. This is consistent with 
the data. The new particles affect the decay widths of h° only at the loop level. Hence 
they can have significant effects only on the loop induced Higgs decays like h° — > 77 and 
hP — > which are not present at the tree-level in the MSSM. 

At the time of the new boson discovery announcement, both ATLAS and CMS exper- 
iments at the LHC reported an excess of events in the pp — > h° — > 77 channel with respect 
to the SM expectations [38]. After including the full 2012 dataset into the analysis, the 
ATLAS collaboration continues to observe an increased signal strength of 1.65^Q'g in the 
diphoton channel [40], while the signal strength measured by the CMS collaboration has 
decreased to 0.78j^2g or Lll^g 2 ,, depending on the analysis method [39]. Naively averag- 
ing the ATLAS and CMS results, we obtain a signal strength of 1.14 ± 0.21 or 1.37 ± 0.22, 
depending on which of the CMS results is used. 

In the following we will discuss the effects of the new leptons and sleptons on the 
h — )• 77 decay and on the mass of the Higgs. We focus on conditions to obtain a moderate 
enhancement of the diphoton rate, which is well in agreement with the data. 
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4.1 h° -»• 77 Width 

The matrix of the leptons (sleptons) to hP couplings is given by the gradients of the lepton 
(slepton) mass matrix with respect to the Higgs VEVs, projected onto the h° eigenstate: 



Yi 



Eh° 



dM% dM% 

sin/3— -^+cos/3— -M, 

v ■ r 9Me j. r 9Me 



dv u 



dv d 



(4.3) 
(4.4) 



where M E is the slepton mass matrix (2.3) while Me is any one of the off-diagonal 
2x2 blocks of the lepton mass matrix (2.2). We also have used the decoupling relation 
a = (3 — tt/2. Therefore, the couplings of h° to the new charged slepton and lepton mass 
eigenstates are given as 



C Eh o = UlY M U 8 , 
C Eh o = uIYeUr , 



(4.5) 
(4.6) 



where Us is the unitary matrix which diagonalizes M E , while Ul and Ur are the unitary 
matrices that diagonalize M.eM e and M. e Me, respectively. 

The new charged leptons and sleptons, Ei and Ei, contribute to the h decay to two 
photons at the one loop level. In the mass basis, the contributions to the decay are given 
by 



h— >77 



X 



Ai(r w ) + -A 1/2 (r t ) + ^ -j^A^te,) + - g ^ ._ 



C Eh° V . , . 



t=l,2 



(4.7) 



where the loop functions for spin 0, spin 1/2 and spin 1 particles are given by [41] 

f(r)-r 



Mr) 

Al/2(T) 

Mr) 



2(r+(r-l)/(r)) 



2 3 3(2r-l)/(r) 



for spin . 
for spin - . 
for spin 1 . 



with 



m h 
4ml 



f(r) 



arcsm 



-iir + log 



1 + Vl 



— l 



1-vT 



for r < 1 . 
for r > 1 . 



(4.8) 
(4.9) 
(4.10) 

(4.11) 
(4.12) 



and m x is the mass of the particle running in the loop. It is instructive to consider the 
asymptotic values of the loop functions for r x <C 1, i.e. when the new particle masses are 
much larger than the half of the lightest Higgs boson mass. Asymptotically A]/ 2 (t — > 0) = 
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4/3 and A (t -> 0) = 1/3, while A 1/2 (r) > 4/3 and A (t) > 1/3 for < r < 1. Note in 
particular that the SM contribution for to/j = 125 GeV is A\{t w ) = —8.3 from the W-boson 
loop and ^Ax/^ijt) = 1.8 from the top quark loop. 

Since the new leptons don't affect the Higgs production channels, the effect on the 
di-photon search channel at the LHC is fully described by the ratio. 

R a(pp -> h) Tjh -» 77) = T(h -> 77) 
77 &sm(pp -> fe) r(/i -> 77)sm r(/i-^77)sM - 

In the limit of rri£ , m e , ^ — )• and vanishing soft SUSY-breaking mass terms, the pre 
factors C Eh o v/Mi and C Eh ov/M^. in (4.7) go to tan/3 and cot/3, respectively, which are 
positive numbers. This leads to destructive interference between the dominant W boson 
contribution and the charged lepton and slepton loops, thus i? 77 < 1 is expected in this 
limit. 

In order to understand how i? 77 > 1 can be achieved, we note that in the asymptotic 
limit the charged slepton contribution to the amplitude can be written as 

vCb, 1 // . „ d \ / „d 



A, = ^ (0)^^4,((sin/ 3 A logdet ^|) + ( cos ^ logdet ^|)) 



(4.14) 



Similarly for leptons, 



A E = A 1/2 {0)Y,^ hlL = ^v((sm/3-^\ogdetME\ + ( cos /^ lo S detM E \\ ■ 

(4.15) 

The second equality is a consequence of the Higgs low energy theorem [42-45], but can 
also be understood by noting that the trace of a matrix is basis independent and using 
log det Ai = tr log Ai . It is useful to evaluate the derivative partially, which leads to 



v 



d det Ml d det Ml 



Ap= — " , I sin/3 ^+cos(3 ^ , (4.16) 

3(detM%) \ dv u dv d / ' v 1 



AE= w^>r fi ^^ +co ^^^)- <4 - 17) 

To obtain an enhancement of i? 77 , it is clear A E + A E must be negative in order to have 
constructive interference with W boson loop. Furthermore it is evident that a large mass 
for all new particles will lead to a suppression of the effects due to the determinant in 
the denominator. This can be used to study the effects of sleptons and leptons separately, 
since either sector can be decoupled by introducing a sufficient amount of supersymmetry 
breaking. 

Conditions for constructive interference from vector-like leptons were studied in detail 
in [1]. A new effect that enters in the case of supersymmetry is the dependence on tan/3, 
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since opposite chirality leptons can't couple to the same Higgs doublet. It follows that 

A E = 4 - , ( 4.i 8) 

3 (1 + t a n^)(-M L M E + tan ^f ) 

Both tan/3 = and tan/3 = oo lead to a suppression of A E , while the effect is maximal 
for tan /3 or order one. Even then, the effect on i? 77 is reduced compared to a non- 
supersymmetric model where leptons of both chiralities are allowed to couple to the same 
Higgs. 

The slepton effects are slightly more involved. Both the vector mass terms Ml, M e 
and the \i term appear in off-diagonal elements of the mass matrix, such that in general 
all four charged slepton fields will mix heavily. 

Absence of tachyons implies that det M% must be positive, therefore the derivatives 
in (4.16) must be negative in order to obtain an enhanced h — > 77 rate. In order to see 
how the different terms in the mass matrices influence As, it is instructive to consider the 
large tan/3 limit. Then all terms proportional to Vd can be neglected. Further setting the 
soft breaking parameters to zero, one obtains 

a 1 2v 2 u y' 2 (M 2 M 2 + 2(M| + jgjigjg + Svjy'^ 2 , . 

E 3 M%Mi-vly>^Ml + vlyf){Ml + vly»V K ' ) 

<h Ml. (4 20) 

3 M 2 M 2 - v 2 y' 2 ^ 2 ' 1 ' ' 

In the last line we have further set y" c — > 0. This result highlights the importance of /j, 2 for 
obtaining a large contribution to i? 77 . While the numerator grows linearly with fj?, the 
denominator can be held roughly constant by appropriately adjusting M E and M E , such 
that the slepton masses are in agreement with the LEP limits. Very large values for Ajj 
can therefore be obtained at the expense of some tuning between the mass parameters. 

In section 5 we will see that vacuum stability considerations will lead to an upper 
bound on fi 2 in the presence of large Yukawa couplings y' c and y" c and small slepton masses, 
which will limit the amount by which i? 77 can be enhanced. Numerical results that cover 
the full parameter region of our model and show the maximal possible enhancement are 
shown in section 6. 

4.2 One Loop Corrections to the h° mass 

To compute the one loop corrections to the h° mass, we use one loop effective potential 
approximation. A supermultiplet i with scalar of mass Ms^ and fermion of mass M Ei 
contributes 




(4.21) 
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Figure 3. Contributions to the Higgs mass Am h o in the m so f t — /i plane for tan/3 = 2 (left) 
and tan/3 = 60 (right). The red (solid) curves are contours of constant Am h o while the green 
(dashed) contours show the mass of the lightest slepton state. The remaining parameters are 
M L = M E 
sleptons. 



200 GeV, y' c — 0.9, y" = 0.7, and we assume equal soft masses mf = m^ oft for all 



to the one-loop effective potential. Here Q is the renormalization scale. Then, as given in 
[30], the correction to the tree- level m? is 



Am 2 



h° 



/sin 2 /? 



dvl 



id 

v u dv u 



+ 



cos 2 /3 



d 2 
dv 2 



1 d 

Vd dv d 



+ sin j3 cos j3 



d 2 

dv u dv d 



i 

(4.22) 



A notable modification of i? 77 requires at least some of the new leptons or sleptons to 
be light. In the absence of soft breaking terms, only fi will induce SUSY breaking in the 
new lepton sector, and therefore corrections to the Higgs mass will remain small. Larger 
corrections are possible if soft terms induce a sizable mass splitting between the sleptons 
and leptons, and it is worth asking whether these corrections can improve upon the fine 
tuning in the MSSM. 

Compared to the corrections typically obtained from top-stop splitting, we can expect 
more moderate contributions due to the absence of a color factor and because the Yukawa 
couplings, that enter in the fourth order in the Higgs mass, are smaller than the top Yukawa 
coupling. In Fig. 3 we show contours of Am h o in the m so ft — A* plane. For these plots, we 
have assumed a base Higgs mass of niQ = 120 GeV in the MSSM and we show the correction 
+ Am^ — mo, with Am 2 ^ coming from (4.22). The sleptons are taken to 



Am h o 



have equal soft masses m 



m 



soft 



, while the soft breaking a and b terms are set to zero. 
The fermionic mass scale is set to M L = Me = 200 GeV and we use y' c = 0.9 and = 0-7, 
close to their fixed point values. 

It is easy to see that even for TeV scale slepton masses, the Higgs mass is lifted by 
at most one GeV, while the corrections are negligible when both leptons and sleptons are 
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close to the weak scale. It is possible to obtain larger corrections by adding large a terms to 
the SUSY breaking Lagrangian. However these soft terms can potentially destabilize the 
vacuum beyond what will be discussed in the next section, so we refrain from introducing 
them. 

Clearly the new lepton sector can not alleviate the fine tuning problem in the MSSM, 
and furthermore will not be sufficient to lift the Higgs mass to 125 GeV in the low tan/3 
case. Instead the necessary contributions could arise from SUSY breaking in the vector-like 
quark sector that is not being discussed in this paper (see e.g. [46] for a recent discussion). 



5 Vacuum Stability 

To analyze the vacuum structure of our model, we need the full scalar potential for the 
Higgs scalars and the new sleptons. We assume that all other scalar fields have masses 
large enough to avoid any stability problems. The scalar potential is obtained from the 
superpotential as 

V = Y1 +l^29 2 aC^2(PiT a <pi) 2 + V soit + radiative corrections. (5.1) 

<t>i a <t>i 

where the sum runs over fa £ {h u , h^, e' L , e' R , e" L , e R }, while fields that are not in danger 
of acquiring a VEV, like the MSSM superpartners or the sneutrinos, are omitted. For the 
different contributions to the potential, we obtain 

Vf = \y' c e' L ~e' R - »hl\ 2 + \y'^el - »h° d \ 2 + \M L E R + y' c h° d ~e' R \ 2 (5.2) 
+ \M E ~e'i + y' c e' L h d \ 2 + \M L e' L + y'Xe" L \ 2 + \M E E' R + y'g'^f , 

Vn = \ 9 l {{hlf - (fc°) 2 + e'l - l" 2 ) 2 + \gl («) 2 - (h° d f - e'? + ~e' R 2 + 2lg - 2~eff , 

(5.3) 

Vr, A = \{gl+gl)&H{hi)\ (5.4) 
^soft = m 2 Hu (h° u ) 2 + m 2 Hd (h° d ) 2 + B»h° u h d + U sofM . (5.5) 

For simplicity we have assumed that all parameters and all fields are real. The Higgs soft 
mass parameters rn 2 H ^ and "rnjj d are usually replaced by the parameters v and tan (3 that 
characterize the electroweak symmetry breaking minimum. For large tan f3, the correction 
Sh must be roughly 5h ~ 1 in order to obtain the correct Higgs mass for the SM-like Higgs 
in the decoupling limit, while for tan f3 ~ 2 a value of 5h ~ 2.5 is necessary. 

As long as the slepton mass matrix has only positive eigenvalues, the above potential 
will have a local minimum characterized by v = 174 GeV and tan /3, with a value for the 
potential V(v , tan f3) = Vq. Additional minima with non-zero VEVs for some of the charged 
slepton fields can be induced by the trilinear terms in (5.2) [47]. 

When the sleptons are heavier than the electroweak scale v = 174 GeV, these addi- 
tional minima typically have a potential energy larger than the electroweak minimum, and 
therefore do not pose a problem. However when the slepton masses are of order v or below, 
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Figure 4. Vacuum stability constraints on the fi parameter. Left: Allowed values for /i in the 
M — y' c plane in the two slepton approximation (5.7). Right: Constraints in the ji — M plane 
for y' c — 0.9, y" = 0.7 and tan/3 = 60. Regions below the solid blue (red) lines are consistent 
with absolute vacuum stability with four (two) light sleptons. The dashed red line indicates the 
approximation (5.7), while the orange region is excluded by the LEP limits on lepton and slepton 
masses. 



some of the charge breaking minima might be lower than the electroweak minimum. Since 
relatively light sleptons are required to obtain a large enhancement of i? 77 , it is clear that 
the vacuum structure of the model must be analyzed carefully. 

Let us first consider the absolute stability condition, namely the condition that (v, tan j3) 
is the global minimum of the potential. In order to see the effects of the various parame- 
ters on the vacuum stability, it is instructive to derive an analytical result in the limit of 
tan/3 — > oo (vd — > 0) and e^ 2 , e"^ — > 0, corresponding to a scenario where large soft mass 
terms for the mirror sleptons prevent them from acquiring a VEV. Further neglecting soft 
a and b terms, the following term is added to the MSSM Higgs potential: 

V = {y'Ae'R? - W u y' c e' L e'R + M 2 L e'l + M 2 g£ + D-terms . (5.6) 

The instability is induced by the trilinear term proportional to y' c fi, whereas all other 
terms are strictly positive. A necessary (but not sufficient) condition for a deeper charge 
breaking minimum is V < 0. Therefore V > is sufficient to guarantee stability of the 
Higgs potential in this limit. After some manipulations, we find 

(5.7) 

where k = — j-^. If the Yukawa coupling is close to the RGE fixed point, this bound 
constrains fi to be at most of the same order as the vector masses. Allowed values for [i 
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as a function of the Yukawa coupling y' c and of the vector mass scale M = Ml = Me 
are shown in Fig. 4 (left). An sizable enhancement of the di- photon rate from the fermion 
sector requires y' c ~ 0.9 and M ~ 200 GeV, which roughly translates to fj, < 250 GeV. 

When all slepton fields are included, analytic bounds on the model parameters become 
impossible to derive. The absolute stability condition can however easily be implemented 
using numerical minimization routines. An additional constraint on \x comes from the LEP 
limit on charged sleptons. 

The right panel of Fig. 4 shows the different constraints in the fi — M plane, for 
tan/3 = 60, y' c = 0.9, y" = 0.7 and Ml = Me = M. The lighter (darker) orange region is 
excluded by the conservative (optimistic) LEP limits. Regions below the solid red line are 
consistent with absolute stability when only two of the sleptons are light, and we see that 
the analytic approximation (dashed) agrees well with the numerical solution. When all 
four slepton fields are light, the stability constraint on fi becomes stronger, as indicated by 
the solid blue line. Overall it is clear that both the LEP limits and the stability constraints 
have to be taken into account when values of /x ~ M are being considered. 

Since the stability limit on \i is more constraining than the LEP bound in many cases, 
it is worth noting that absolute stability is not a necessary condition for the model to be 
viable phenomenologically. After all, the measured Higgs mass of 125 GeV implies that the 
standard model itself is only meta-stable [48, 49]. This means that while the electroweak 
minimum is not the global minimum of the radiatively corrected Higgs potential, the tun- 
neling rate to the true global minimum is suppressed enough to guarantee a lifetime larger 
than the age of the universe. 

Imposing the weaker meta-stability bound on our model could allow larger values of 
\i which in turn will lead to higher attainable diphoton rates, as was discussed in Sec. 4. 
If the electroweak minimum is not the global minimum of the potential, the probability to 
transition into the true charge breaking vacuum per unit volume and time depends on the 
decay rate [50] 



Here A is a dimensionful parameter that is expected to be of fourth order in the electroweak 
scale, A ~ v 4 , and Se is the Euclidean action of the bounce solution corresponding to the 
transition from Higgs vacuum to the charge symmetry breaking true vacuum. The age 
of the universe is about 1.37 x 10 10 years, which implies that Se > 400 is a necessary 
condition for our existence. 

The metastability bound computation is performed using the package CosmoTransi- 
tions [51]. It finds the bounce solution for the transition from the false vacuum to the real 
vacuum of a multi-dimensional scalar potential by the method of path deformation. The 
results of this computation and its impact on h — > 77 enhancement are presented in Sec. 6. 
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6 Results 



The model has several distinct regions of parameter space that can lead to an enhanced 
i? 77 . In the following we will discuss the most interesting scenarios and their strengths and 
weaknesses. 

6.1 Vector-like leptons 

First we will consider the scenario where only the new leptons are light, while the sleptons 
are lifted to the TeV scale by the diagonal soft mass terms in (2.3). TeV scale superpartners 
are sufficient to protect the Higgs quartic coupling from running to negative values below 
the scale where supersymmetry is restored. Furthermore large mass terms for the sleptons 
protect this scenario from vacuum instabilities and ensure that the electroweak symmetry 
breaking vacuum is a global minimum. 

The conditions for obtaining a non-negligible, positive contribution to i? 77 are sum- 
marized in Eqn. (4.17). First, we require that the combination M^Me — v 2 y' c y" c sin j3 cos /3 
is positive but not too large, in order to obtain the correct sign for Ae- Furthermore 
tan/3 must be of order one, otherwise the contribution will be suppressed by l/tan/3. The 
latter condition emerges because the effective Yukawa couplings of both leptons and mirror 
leptons to the lightest CP even Higgs boson are rescaled by sin /3 and cos (3 respectively, so 
their ratio should not be too small. 

One immediate concern is that a mass of 125 GeV for the lightest Higgs boson is 
difficult to obtain with tan/3 of order one within the MSSM. While we have seen in Sec. 4 
that the new slepton sector does not significantly improve the situation, there are other 
ways to lift the Higgs mass without affecting it's low energy phenomenology. The most 
straightforward solution is to assume that additional one loop contributions beyond those 
of the top quark come from the vector-like quark sector that accompanies the leptons if 
they are implemented in complete SO(10) multiplets. Other possibilities include scalar 
singlet extensions of the MSSM or additional gauge interactions [27]. For the remainder of 
this section, we will assume that one of these mechanisms is at work, but does not otherwise 
affect the phenomenology of the lightest Higgs boson. 

Compared with the results of [1], we expect that the enhancement of i? 77 is suppressed 
roughly by l/tan/3. Note that tan/3 > 1.5 is required to ensure perturbativity of the top 
Yukawa coupling in the presence of sizable Yukawa couplings for the new leptons. In Fig. 5 
left we show the enhancement that can be obtained for tan/3 = 2 and y' c = 0.7, y" = 0.9, 
close to their respective fixed point values, in the Ml — Me plane. While the features of 
the plot are similar, the maximal enhancement that can be obtained for a lightest lepton 
mass of order 100 GeV is i? 77 < 1.3, compared with i? 77 < 1.6 in the non-supersymmetric 
case. 

To further illustrate the importance of tan /3, in the right panel of Fig. 5, we show the 
contours of i? 77 in the plane of tan /3 and m ei , the mass of the lightest new lepton. For 
masses around 100 GeV the enhancement is reduced from around 30% at tan j3 = 2 to 
below 10% for tan /3 = 10. 
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Figure 5. Enhancement of the h — > 77 rate from vector lcptons only, with charged Yukawa 
couplings Yj. — 0.9 and Y" c " = 0.7. Left: Contours of R 11 in the Ml — Me plane (blue, dashed), 
for tan/3 = 2. Also shown are contours for the mass of the lightest charged lepton state (orange, 
solid). Right: Enhancement i? 77 as function of the lightest charged lepton mass for different 
values of tan j3. 

Overall one can see that the cost of imposing GUT scale stability cuts the enhancement 
of i? 77 in half. If these new fermions are the only particles that have significants effects on 
the Higgs phenomenology, only a modest enhancement of i? 77 of around 30% is compatible 
with vacuum stability, perturbativity, and grand unification. 

6.2 Supersymmetric leptons 

Here, we will first assume that at the tree level the only source of SUSY breaking in the 
new lepton sector is through the /i term, such that a minimal number of parameters are 
added to the MSSM. When \i = each lepton is accompanied by two sleptons with the 
same mass. These degenerate states are split when jx is nonzero. Therefore, in the absence 
of soft breaking parameters, the mass of the lightest slepton is always lower than the 
corresponding lepton mass. 

Let us first consider the large tan j3 limit, where only the sleptons contribute signifi- 
cantly to Ryj. As before, the Yukawa couplings are taken to be y' c = 0.9 and y" = 0.7, 
and we will usually take Ml = Me = M. Since all sleptons are now light, in addition to 
the LEP bounds we have to impose vacuum (meta-) stability for the model to be viable. 
Fig. 6 shows the diphoton rate enhancement that can be obtained for tan/3 = 60. Stable, 
meta-stable and unstable regions of parameter space are indicated in white, yellow, and 
red, respectively. The orange shaded region is excluded since we required m^, > m/j/2. 
Contours of constant it! 77 are shown (blue, dashed) as well as contours of constant 
(green, solid). The left plot, where \x = 280 GeV, demonstrates the usual result that 
M L = Me maximizes the enhancement, and we can see that an enhancement of ii 77 by 
40% is possible for ms > 90 GeV and meta-stability. It is evident from the structure 
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Figure 6. Diphoton enhancement rate with light leptons and sleptons, for tan/3 = 60. Regions 
with stable, meta-stable or unstable vacuum are shaded white, yellow and red, respectively. Regions 
excluded by LEP are shaded orange. The blue (dashed) lines are contours of constant i? 77 . The 
brown (solid) lines are the slepton mass contours of 90, 150 GeV for the left and 62.5, 90, 150 GeV 
for the right panel reading away from the unstable region. The left panel uses \i = 280 GeV. 

of the lepton mass matrix given in Eq. (2.2) that the lepton masses are well above the 
conservative LEP bound of 100 GeV for the region of the parameter space shown here. 

The right plot of Fig. 6 shows the interplay of the different constraints as \x is increased. 
Initially, the LEP bound is more constraining, and enhancements of at most 40% are 
possible with > 90 GeV. Following the m^, =90 GeV contour to higher values of \i 
we note that i? 77 increases up to 1.5 before we eventually hit the stability bound. Lowering 
the lightest slepton mass is only possible for /i < 280 GeV. In that region, enhancements 
of 100% or more seem possible when m& is very close to the absolute lower bound of 
62.5 GeV. 

In Fig. 7 we show the diphoton rates for tan /3 = 2, while all other parameters are the 
same as before. Here we expect that the fermions contribute up to 30% enhancement to 
the diphoton rate, such that higher total rates should be possible. However it turns out 
that the vacuum stability constraint forces us to consider lower values of fj,, such that the 
overall enhancement is not significantly higher than in the large tan f3 case. 

More precisely, it can be seen in the left panel that an enhancement of over 50% is 
possible with \x as low as 175 GeV and m_g > 90 GeV, which is not possible in the large 
tan/3 case. However the right plot clearly shows that larger values of n do not further 
increase i? 77 , since the stability constraint requires larger values of M at the same time. 
Moving along the mj, =90 GeV contour, the ratio of fermionic to scalar contributions, 
Ae/AE, decreases from about 0.5 at \i = 120 GeV to 0.3 at \x = 175 GeV. Enhancements 
of order 100% are again only possible if we lower the slepton masses below the conservative 
LEP bound. 
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Figure 7. Same as Fig. 6 for tan/3 = 2. In the left panel fj, = 175 GeV. 



Since the soft mass terms are set to zero, the vector masses Ml and Me are needed 
to lift the sleptons above the LEP bound. We have seen in Sec. 6.1 that the leptonic 
contribution is maximized around Ml = Me = 200 GeV. However such low values are 
not allowed since the mass of the lightest slepton would drop below the LEP limit. It is 
therefore interesting to explore what happens when we add small soft terms to stabilize 
the slepton masses, such that both the leptonic and the slepton contributions to i? 77 can 
be maximized. 

It is instructive to qualitatively discuss the impact of the various elements of the slepton 
mass matrix before proceeding to finding the maximum possible i? 77 in the presence of soft 
terms. The main role of the off-diagonal terms is to increase the difference between the 
eigenvalues, while the sum of the eigenvalues is given by the trace of the mass matrix. 
Therefore, the parameters that only appear in off-diagonal terms increase the difference 
between eigenvalues keeping the sum constant, effectively resulting in lowering the lowest 
eigenvalue. 

It is clear from the structure of the slepton mass matrix given by Eqn. (2.3) that the 
fj, parameter and the bilinear holomorphic soft terms are only present in the off-diagonal 
terms. Thus both of these can drive the lightest slepton mass below the LEP bound as 
well as destabilize the vacuum. On the other hand, as illustrated in Eqn. (4.19), these 
off-diagonal terms are essential to produce i? 77 significantly higher than 1. These two 
parameters have similar impact on i? 77 , lightest slepton mass and stability. Therefore, in 
order to find the maximum enhancement, we set the holomorphic soft terms to zero and 
achieve the h — > 77 enhancement only through /x as before. 

The non-holomorphic soft terms mf that appear on the diagonal of (2.3) can then be 
used to lift the lightest slepton mass above the LEP limit. For simplicity we will take 
them to be equal, mf = wnt & . Fig. 8 shows the enhancement rates that can be obtained 
for tan/3 = 60 (left) and tan/3 = 2 (right). The vector mass scale M is chosen such that 
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Figure 8. Diphoton rates with nonzero soft mass terms for the sleptons, mf = m^ oft . The vector 
mass terms Ml = Me = M are chosen such that the lightest lepton has a mass around 100 GeV. 
The left panel shows the result for tan/3 = 60 and M =150 GeV, while the right panel is for 
tan (3 = 2 and M = 190 GeV. Colors same as in Fi.g 6. 



the lightest lepton mass is about 100 GeV, corresponding to M = 150 GeV (tan/3 = 2) 
and M = 190 GeV (tan/3 = 60), respectively. In the case of large tan/3, we note that an 
enhancement of 50% is now possible without going below the conservative LEP limit on the 
slepton masses, an improvement of about 10% compared to the case without soft terms. 
Similarly, in the low tan/3 case, we can now get up to 75% enhancement of the diphoton 
rate without making the sleptons dangerously light. Furthermore one should note that 
the absolute stability limit is less constraining now, and values of i? 77 of 1.5 and 1.6 are 
compatible with absolute stability for high and low tan /3, respectively. 

6.3 Split sleptons 

In Sec. 5 we have seen that imposing absolute vacuum stability puts modest constraints 
on the parameter space when only two charged sleptons are allowed to get a VeV, but that 
the constraints become strong when all four sleptons are taken into consideration. 

Here we will consider a scenario where the sleptons are split by TeV scale soft masses 
m|. for the double primed fields e" L and e" R . This will improve the vacuum stability, since 
non-zero expectation values for these fields are unlikely to give vacua deeper than the 
EWSB vacuum. The remaining light charged degrees of freedom are the two sleptons e' L 
and e! L as well as both the charged leptons. Phenomenologically this is a combination of 
the vector- like lepton model [1] with the light stau scenario [34]. 

Absolute vacuum stability is guaranteed provided that the inequality (5.7) is satisfied. 
As can be seen from Fig. 4 the conservative LEP limit on the mass of the lightest charged 
particle, m_g > 90 GeV, is in general more constraining, so that the vacuum stability 
constraint is automatically satisfied for most phenomenologically viable parameter points. 
Regions of parameter space that are not absolutely stable will be indicated in the plots, 
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Figure 9. Enhancement of the h — > 77 rate in the split slepton scenario for tan/3 = 2 (left) and 
tan/3 = 30 (right), with charged Yukawa couplings Y' c — 0.9 and Y" — 0.7. Shown are iso-contours 
of i? 77 (blue, dashed), of the lightest slepton mass (green, solid) and of the lightest lepton mass 
(red, dotted). Note that values of /x below 100 GeV are in conflict with limits on chargino masses. 
The parameter space outside of the absolute stability region is shaded yellow. We do not show the 
meta-stability region since it lies outside of the LEP allowed region. 

but one should keep in mind that they can still be phenomenologically viable if the much 
weaker meta-stability bound is satisfied. 

We can again distinguish two scenarios, low tan j3, where leptons and the light sleptons 
contribute to i? 77 , and the large tan/3 regime, where only the two light sleptons can make 
a notable contribution. 

Let us first consider the small tan f3 case. For defmiteness we take tan j3 = 2 and, as in 
the previous sections, y' c = 0.9 and = 0.7. In Fig. 9 (left) we show i? 77 in the fi-M plane, 
where M = Ml = Me is the vector mass scale. As can be seen, the enhancement can 
be increased both by lowering the lightest slepton mass (i.e. increasing fi) or by lowering 
the lightest lepton mass (decreasing M). When both masses are close to 100 GeV an 
enhancement of about 40% can be obtained. Values up to 70% can be reached by lowering 
the lightest slepton mass further, while still being consistent with limits from the LEP 
experiments and with absolute vacuum stability. 

When tan /3 is increased, the leptonic contributions are suppressed. The case of tan /3 = 
30 is shown in Fig. 9 (right). i? 77 is now independent of the lightest lepton mass, and an 
enhancement of 30% or more can only be obtained when the lightest slepton mass is below 
100 GeV, and only for sufficiently large values of fi. Imposing absolute vacuum stability is 
more constraining here, which disfavors the large fi region where i2 77 can be of order 1.5 
or higher. 

In Sees. 6.1-6.3, we have focussed on the maximal possible values of i? 77 that can be 
obtained in each scenario. In contrast to many existing models, we find that at least in some 
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scenarios it is possible to obtain enhancements of more than 50% while at the same time 
the new particle masses can be kept above the conservative LEP bound, vacuum stability 
is maintained and all couplings remain perturbative up to high scales. Nevertheless it is 
evident from our figures that some amount of tuning of the parameters is necessary to obtain 
such large values for i? 77 , while more modest enhancements of 20% — 40% are possible in 
larger regions of parameter space and thus appear more natural. Such an enhancement is 
well in agreement with the signal strength indicated by the combination of the updated 
ATLAS and CMS results [39, 40]. 

7 Conclusions 

The recent discovery of a Higgs-like particle at the LHC opens a new era in particle physics. 
A very important task is to study the properties of this particle in detail, and to analyze 
any possible deviation from the SM predictions that might signal the presence of new 
physics. Currently, the measured Higgs-induced diphoton production rate is 2.3 a above 
the SM prediction at the ATLAS experiment [40]. This provides a motivation for the study 
of new physics scenarios that can lead to an enhancement of the diphoton decay width. 
Vector- like leptons provide an extension of the SM that leads to such an enhancement. 

Quite generally, the presence of new weakly interacting particles with strong couplings 
to the Higgs boson can provide an enhancement of the loop-induced diphoton rate, but 
also lead to the presence of new vacua deeper than the physical one. In the case of vector- 
like leptons such vacua arise at large values of the Higgs fields due to the evolution of the 
quartic coupling of the Higgs to negative values. Enhancements of the diphoton rate to 
values larger than 1.5 times the SM one can only be obtained if new physics stabilizes the 
vacuum at scales smaller than a few TeV. Supersymmetry provides a natural extension of 
this model in which the vacuum stability problem of the Higgs sector is provided by the 
contributions of sleptons. 

In this article we study the supersymmetric extension of the vector-like lepton theory 
introduced in Ref. [1]. We showed the inclusion of a whole vector-like family can lead to 
a unified theory with values of the gauge couplings that are close to the non-perturbative 
bound at scales close to the GUT scale. In order to enhance the effects on the Higgs 
diphoton decay rate, we chose values of the Yukawa couplings leading to large values at 
the GUT scale, but still consistent with a perturbative treatment of the theory. Vector-like 
squarks are considered to be heavy and therefore have an impact only in the determination 
of the Higgs mass. 

The phenomenological properties of this supersymmetric extension depends strongly on 
the values of the soft breaking parameters. For large values of the scalar soft supersymmetry 
breaking parameters, the theory at low energies reduced to the one studied before. However, 
for the same values of the Yukawa couplings, the lepton contributions are suppressed by a 
sin 2/3/2 factor and, together with modified values of the perturbativity bounds on these 
couplings with respect to the SM ones, the enhancements of the diphoton rate remain 
smaller than about 30 percent for lepton masses above 100 GeV. 
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For small values of the scalar soft supersymmetry-breaking parameters, the main source 
of supersymmetry breaking in the Higgs-slepton potential is provided by the Higgsino mass 
parameter fj,. For light sleptons, large values of fi tend to induce new charge breaking 
minima in the spectrum, and therefore in the region consistent with vacuum stability light 
sleptons are associated with relatively light leptons. It follows that both the fermion and 
the scalar lepton contributions to the Higgs-induced diphoton production rate tend to be 
important, except for the large tan/3 regime where the lepton contributions decouple. We 
find that for lepton and slepton masses larger than 100 GeV enhancements of order 40 
and 50 percent may be obtained for small and larger values of tan/3 (tan/3 = 2 and 60), 
respectively. More generally, since the LEP bound depends on the value of the neutrino 
mass parameters, one can consider leptons and sleptons as light as 62.5 GeV, for which 
much larger enhancements of the Higgs diphoton rate may be obtained. 

Finally, we considered a split slepton scenario, in which the soft supersymmetry- 
breaking parameters of the new right-handed leptons are considered to be large, while 
the left-handed ones are kept small. In such a case, the theory is similar to the light-stau 
scenario, but the lepton contributions remain relevant. We showed that in such a case 
enhancements of the Higgs-induced diphoton rate of the order of 50 percent and 30 percent 
can be obtained for small and large values of tan /3 for lepton and slepton masses above 
100 GeV, while as before larger values may be obtained if these bounds were relaxed. 

In order to avoid flavor problems, we have introduced a parity symmetry under which 
the new states are odd. If in addition R-parity is conserved this guarantees three stable 
particles: the ordinary MSSM LSP as well as the lightest parity odd leptons and slep- 
tons. Assuming that in each sector a neutral particle is the lightest state, this leads to 
a scenario with multicomponent dark matter with possibly interesting phenomenological 
consequences. 

Even in the light of the recently presented CMS results [39], there are large regions 
of parameter space in which the vector like leptons and sleptons are present at the weak 
scale and remain compatible with data. The masses of the new particles in such region can 
be small enough to make them interesting to study from the perspective of dark matter 
and low-scale leptogenesis. They can also be probed at the LHC, what serves as another 
motivation to study supersymmetric vector-like leptons in detail. 
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Note added: While finalizing the manuscript for submission, Ref. [16] appeared, which 
also considers a supersymmetric extension for the vector-like lepton scenario. Different 
from the present work, the authors do not impose perturbativity of the Yukawa couplings 
at the GUT scale, and therefore allow larger Yukawa couplings at the electroweak scale. We 
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also have shown that the relevance of vacuum stability constraints depends on the choice 
of soft breaking parameters. 
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